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The short and intermediate range order of an amorphous GeSe4 alloy produced by Mechanical 
Alloying were studied by Reverse Monte Carlo simulations of its x-ray total structure factor and 
Raman scattering. The simulations were used to compute the G™Gc( r )> C™gj,(r) and Gg^j(r) 
partial distribution functions and the Sq^}q c (K), Sq™c? c (K) and Sg^£.(K) partial structure factors. 
We calculated the coordination numbers and interatomic distances for the first and second neighbors 
and the bond-angle distribution functions O i:) ;(cos 9). The data obtained indicate that the structure 
of the alloy has important differences when compared to alloys prepared by other techniques. There 
are a high number of Se-Se pairs in the first shell, and some of the tetrahedral units formed seemed 
to be connected by Se-Se bridges. 

PACS numbers: 61.10.Eq, 61.43.Bn, 05.10.Ln, 87.64.Je 



I. INTRODUCTION 



Amorphous semiconductor materials like chalcogenide glasses present a great potential for application in techno- 
logical devices, such as optical fibers, memory materials and switching devices, but their use is limited due to several 
factors. One of them is the difficulty in obtaining information about atomic structures, which define the short-range 
order (SRO) of the alloy. In this context, the structures of amorphous Ge^Sei-a; (a-Ge^Sei-^) and liquid Ge^Sei-^ 
(Z-Ge^Sei-a:), in particular Ge33Se67 (GeSe2), have been extensively studied by several experimental techniques, like 
neutron diffraction (ND)i2i2i^, x-ray diffraction (XRD)&£, extended x-ray absorption fine structure (EXAFS)^ and 
Raman spectroscopy (RS)i2*ii*i£. On the theoretical side, molecular dynamics simulations (MD)i2ii2ii£ii£iiLi£ii2i22i2i 
and reverse Monte Carlo (RMC) simulations 7 have been carried out to understand the SRO in these liquids and glasses 
in terms of two possible and distinct models. In the first one the distribution of bonds in the structure is purely ran- 
domic and determined by the local coordination numbers and composition. In the second one, there is a strong SRO 
and the structure is formed by well defined structural units, e.g., corner-sharing GeSe4/ 2 (CS) tetrahedral and edge- 
sharing Ge2Se 8 / 2 (ES) bitetrahedral units. The distribution of these units gives raise to a medium, or intermediate, 
range order (IRO), whose signature is the appearance of a first sharp diffraction peak (FSDP) in the neutroni^i^ 
or x-ray structure factors^* 7 - at many compositions. In particular, ND experiments performed on melt-quenched^*^ 
(MQ) GeSe2 (MQ-GeSe 2 ) showed a FSDP in the total structure factor S{K) which was associated with correlations 
in the range of 5-6 A. As described in Ref. H, this alloy is formed by CS and ES units with heteropolar bonds but 
there are homopolar bonds in very small quantities. It should be noted, however, that almost all available data about 
a-Ge^Sei-a; alloys were determined for MQ samples, and the preparation method can affect the SRO and IRO. We 
have recently verified this assumption in a study about the structure of o-Ge3oSe7oi produced by Mechanical Alloying 
(MA) 22 . In our case, an unexpected large number of Se-Se pairs is found in the first coordination shell, suggesting 
that the tetrahedral units are linked by Se-Se "bridges". Takeuchi et alML, by comparing the structures of a-Ge3oSe7o 
produced by MQ and by vacuum evaporation (VE) techniques, and Tani et al& studying the a-GeSe2 produced by 
Mechanical Grinding (MG) of its crystalline counterpart have also found some structural differences among alloys 
produced by different methods. These differences are important because some physicochemical properties can be 
altered and improved as desired by choosing an appropriate preparation method. 

In this paper, we investigated the SRO and IRO of an amorphous GeSe4 alloy produced by MA (MA-a-GeSe4) 
starting from the elemental powders of Ge and Se using Raman spectroscopy (RS), X-ray diffraction and reverse Monte 
Carlo (RMC) simulations 2 ^* 2 ^* 2 ^* 2 ^ of its XRD S(K). We were interested in studying two main points. First of all we 
would like to know if the alloy produced by MA contains CS or ES units. Besides that, even if these units are formed 
the SRO and the IRO of the alloy can be significantly altered by the high quantity of defects and disorder introduced 
by the MA process when compared to MQ samples, for instance. Therefore, the second point is to determine the 
local structure of the alloy itself, finding coordination numbers and interatomic distances. At our knowledge, this is 
the first time that such study is reported concerning an a-GeSe4 alloy produced by MA. 
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II. THEORETICAL BACKGROUND 
A. Structure Factors 

1. Faber and Ziman structure factors 

According to Faber and ZimanSl, the total structure factor S(K ) is obtained from the scattered intensity per atom 
I a (K) through 

n n 

i=i j=i 

where if is the transferred momentum, Sij(K) are the partial structure factors and Wij(K) are given by 



and 



i 

Here, fi{K) is the atomic scattering factor and Cj is the concentration of atoms of type i. The partial reduced 
distribution functions £?y (r) are related to <Sy (if) through 

2 f°° 

G y ( r ) = ~ K [Sij (K) - 1] sin^r) dif . 
77 Jo 

From the Gy (r) functions the partial radial distribution function RDFy (r) can be calculated by 

RDFjj(r) = AirpoCjr 2 + rG l3 (r) . 

where po is the density of the alloy (in atoms/A 3 ). Interatomic distances are obtained from the maxima of Gij(r) 
and coordination numbers are calculated by integrating the peaks of RDFy (r) . 

2. Bathia and Thornton structure factors 

The Bathia- Thornton (BT) structure factors can be related to the FZ ones2£. For a binary alloy the BT number- 
number structure factor <Snn(-K") is given by 

Ssn(K) = c\S lx {K) + c\S 22 {K) + 2c x c 2 S X2 (K) , (1) 

where Sij(K) are the FZ partial structure factors and Cj is the concentration of element i. The BT number- 
concentration structure factor S^q(K) is 

S NC (K) = c x c 2 {c x [S xx (K) - S X2 (K)] - c 2 [S 22 (K) - S X2 (K)]} , (2) 
and the BT concentration-concentration structure factor Sec (K) is found through 

Scc(K) = Cl c 2 {l + c x c 2 [S xx {K) +S 22 {K) - 2S X2 {K)] } . (3) 
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B. RMC Method 



The basic idea and the algorithm of the standard RMC method are described elsewhereSiS 4 *SSiS& and its application 
to different materials is reported in the literaturo 29 i 30 i 31 i 32 i 33 i 34 i 35 i 36 i 37 i 38 i 39 i 40 i 41 . In the RMC procedure, a three- 
dimensional arrangement of atoms with the same density and chemical composition of the alloy is placed into a 
cell (usually cubic) with periodic boundary conditions and the G^ MC (r) functions corresponding to it are directly 
calculated through 



n RMC( r \ 



4:TTpor 2 Ar 



where n RM (r) is the number of atoms at a distance between r and r + Ar from the central atom, averaged over all 
atoms. By allowing the atoms to move (one at each time) inside the cell, the G R (r) functions can be changed and, 
as a consequence, S^ MC (K) and S RMC (K) are changed. Thus, S RMC (K) is compared to the S(K) factor in order to 
minimize the differences between them. The function to be minimized is 



^ 2 = -^[W)-«s RMC (^)] 2 , ( 4 ) 

i=l 

where the sum is over m experimental points and S is related to the experimental error in S(K). If the movement 
decreases ip 2 , it is always accepted. If it increases %p 2 , it is accepted with a probability given by exp(— At/j 2 /2); 
otherwise it is rejected. As this process is iterated ip 2 decreases until it reaches an equilibrium value. Thus, the 
atomic configuration corresponding to equilibrium should be consistent with the experimental total structure factor 
within the experimental error. By using the G RMC (r) functions the coordination numbers and interatomic distances 
can be calculated. In addition, the bond-angle distributions 0jj/(cos#) can also be determined. 



III. EXPERIMENTAL PROCEDURES 



The MA-a-GeSe4 alloy was produced by considering a binary mixture of high-purity elemental powders of germa- 
nium (Alfa Aesar 99.999% purity, particle size < 150 //m) and selenium (Alfa Aesar 99.999% purity, particle size < 
150 fj,m) that was sealed together with several steel balls into a cylindrical steel vial under an argon atmosphere. The 
ball-to-powder weight ratio was 5:1. A Spex Mixer/Mill model 8000 was used to perform MA at room temperature. 
The mixture was continuously milled for 50 h. A ventilation system was used to keep the vial temperature close to 
room temperature. The XRD pattern was recorded in a powder Siemens diffractometer equipped with a graphite 
monochromator, using the CuK Q line (A = 1.5418 A). The total structure factor S(K) was computed from the XRD 
pattern after corrections for polarization, absorption, and inelastic scattering, following the procedure described by 
Wagner—. The /' and /" values were taken from a table compiled by Sasaki 43 . Raman measurements were performed 
with a T64000 Jobin-Yvon triple monochromator coupled to a cooled CCD detector and a conventional photon count- 
ing system. The 5145 A line of an argon ion laser was used as exciting light, always in backscattering geometry. The 
output power of the laser was kept at about 200 mW to avoid overheating the samples. All Raman measurements 
were performed at room temperature. 



IV. RESULTS AND DISCUSSION 



A. Raman Scattering 



Figure ^ shows the RS spectra of MA-a-GeSe4, c-Se and c-Ge. The alloy has bands at around 195, 215, 237 and 
255 cm -1 . The band at 195 cm -1 is associated with the Ai breathing mode of CS units and the 215 cm -1 band 
is related to the A\ brcathinglike motions (companion peak) of Se in the ES unitsi 1 ^^. The difference in the 
intensities of these peaks indicates that the alloy is formed basically by CS tetrahedra and ES tetrahedra are found 
in a small quantity. The very weak shoulder around 237 cm -1 is associated with A\ and E modes of Se chains, and 
they are also seen in c-Se (see Fig. QJ. The broad and intense band at 255 cm -1 is related to A\ and Ei modes of 
Se„ ringg 8 ' 10 - 4 ^!^ 4 ^. It is important to note that the band at around 165 cm" 1 , which is associated with Ge-Ge pairs 
vibrations in ethanlike units, is not seen in the spectrum of the alloy (as well as the other bands of c-Ge), indicating 
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that the quantity of Ge-Ge pairs is very low in the alloy. In addition, the bands related to Se-Se pairs are not seen in 
the spectra of the alloys produced by MQ 4 ^ or VEi£ techniques at this composition, but they exist in MA-a-Ge 3 oSe7o^ 
and in MA-a-GeSe4, and their intensities suggest that the number of Se-Se pairs in our alloy may be relevant. These 
results indicate that the tetrahedral units are formed during the MA process and there is a preference for CS units. 
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FIG. 1: RS spectra of (a) MA-a-GeSe 4 , (b) c-Se and (c) c-Ge. 



B. X-ray Diffraction and RMC Simulations 



Figure |2] shows the experimental XRD S(K) (full line) for our alloy and the experimental ND S(K) (dashed line) 
given in Ref. for MQ-a-GeSe4. The FSDP is seen at around 1.1 A -1 . It is lower than those shown in Refs. and 
[ij, indicating that the IRO in the alloy produced by MA is less pronounced than in the MQ-GeSe4 samples 1 ' 4 . The 
FSDP is known to be much dependent on Ge-Ge and, to a lesser extent, on Ge-Se correlationsiii 4 * 4 ^^. Therefore 
these correlations have a different behavior in MA-a-GeSe4. This fact is also verified in MA-a-Ge3oSe7(j£. 

S(K) was simulated using the RMC program available on the Internet 25 . To perform the simulations we have 
considered a cubic cell with 16000 atoms (3200 Ge and 12800 Se), S = 0.002, and a mean atomic number density 
po = 0.03834 atoms/A 3 . This value was found from the slope of the straight line (—4irpor) fitting the initial part 
(until the first minimum) of the total G(r) function^. The minimum distances between atoms were fixed at the 
beginning of the simulations at tgc-Gc = 2.0 A, tqc-Sc = 1-8 A and rs c -Sc = 1.7 A. The iS RMC (iQ obtained from the 
simulations (squares) is also shown in Fig. [21 and there is a very good agreement with the experimental S(K). 
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FIG. 2: Experimental (full line) and simulated (squares) total structure factor for MA-a-GeSe4 together with the ND total 
structure factor given in Ref. (dashed line) (for a better comparison it was cut at ii'max = 8.0 A -1 ). 
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FIG. 3: GGe-Ge(r), GGe-Sc(r) and Gsc-Sc(r) functions obtained from RMC simulations of MA-a-GeSe4 (thin solid lines) and 
MA-a-Ge3oSe7o (thick solid lines) and also the ND functions for MQ-GeSe2 (dashed lines, Ref. 0). 

First hard sphere simulations without experimental data were carried out to avoid possible memory effects of the 
initial configurations in the results. Then unconstrained runs (i.e. when only hard sphere diameters and experimental 
data were used during the simulation) were carried out. These runs led to essentially identical G^ MC (r) functions, 
because of the proximity of the atomic numbers and scattering factors of Ge and Se. In the next series of simulations we 
used the coordination numbers we have found for MA-o-Ge3oSe7o^ as starting coordination constraints (A*Ge-Ge = 0.25, 
-^Ge-Sc = 3.5 and As c -Sc — 1-25), which were then allowed to vary freely. The results obtained from the best simulation 
achieved are shown in Fig. [21 and they are discussed below. As a final test, we also tried to make simulations forcing 
the Se-Se coordination to be two, as it is in a-Se, but again the simulations did not show a good convergence. When 
this constraint is released, the Se-Se coordination number decreases until it reaches the number we found for the 
previous case (iVso-Sc = 1.71, see text below). This indicates that there are no amorphous or crystalline selenium in 
the mixture, and the Raman bands at 237 and 255 cm -1 are associated with Se atoms that belong to the MA-a-GeSe4 
alloy. 

1. Pair Distribution Functions 

Figure |3 shows the G^ MC (r) functions obtained from the RMC simulations of MA-a-GeSe4 compared to those 
found experimentally for MQ-GeSe2 using ND with isotopic substitution^. Since the Gij (r) functions shown in Fig. 
3 of Ref. [j for MQ-GeSe4 were calculated and their combination do not reproduce quite well the experimental G(r) 
function for this alloy, we have chosen not to compare our results with those functions. We have also added to Fig. 
the Gy MC (r) functions obtained for MA-a-Ge3oSe7(j£ since it is interesting to compare both alloys produced by MA. 
The intensity of the first peak in Gg^g^ is higher in the present alloy, confirming their existence in a larger quantity 
in the first coordination shell, as indicated by the RS measurement. The two first peaks of the G™Ge( r ) function, 
which correspond to Ge-Ge first and second neighbors, show up around 2.35 and 3.80 A. The first peak occurs at 
a distance a little shorter than that found in MQ-GeSe2^, but the second is displaced towards higher-r values by 
0.23 A. In addition, the peak at around 3.0 A seen in MQ-GeSe2 is not resolved in MA-a-GeSe4. Remembering that 
the distance between two Ge atoms in adjacent ES and CS units are found at 3.02 A and 3.60 A, respectively^^, 
it can be seen that the fraction of ES units in MA-a-GeSe4 is low, again in agreement with the results obtained by 
RS. Since the intensity of the FSDP in S(K) seems to be related to the quantity of ES tctrahedra— the low 
quantity of ES units in our alloy could explain the low intensity of the FSDP. 
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The first peak of Gg^ G (r) function is located at 2.36 A. This shell in our alloy is lower and broader than in 
MQ-GeSe 2 but, due to the difference in densities, Ge-Se coordination numbers are almost the same in both alloys. 
The next peak appears at 3.84 A, and it is higher than that at 4.95 A. In the MQ-GeSe2 samples, there is a peak 
around 3.02 A which is not seen, or at least not resolved, in MA-a-GeSe4, and there are peaks at 3.78 A (smaller) 
and 4.66 A (higher). In c-GeSe2 an ES Ge atom has Se neighbors in the range 4.6 < r < 5.3 A, and a CS Ge atom 
has Se neighbors in the range 4.0 < r < 4.8 A 1 ^ Then, we have associated the peak at 3.80 A with CS units and that 
at 4.95 A with CS and ES units. 

The first peak of Gg^g^(r) function is located at 2.33 A, and it corresponds to a coordination number -/Vso-Sc = 1.71, 
which is much higher than that obtained for the MQ-GeSe2 samples (As c -Sc = 0.20) 3 ^. This suggests that some of 
the tetrahedral units are connected by Se "bridges", forming sequences such as Ge-Se-Se-Se-Ge. As a consequence, 
Ge-Ge pairs should be found at higher distance values. This agrees with the results obtained from the G™Gc( r ) 
function. The next peak appears at 3. 75 A , which gives a ratio of Ge-Se:Se-Se distances of 0.632. The value expected 
for ideal tetrahedral coordination is y/3/8 = 0.612, indicating that the tetrahedral (CS or ES) units are distorted in 
our alloy. It is interesting to note that in MQ-GeSe2 samples no peaks are found from » 4.5 A to « 5.9 A either 
considering ND results 3 ' 5 or MD simulations 14 20 . On the other hand, in MA-a-GeSe4 there are peaks at 4.95 and 
5.75 A, and we believe these peaks are related to the distances between Se atoms in the "bridges" and Se atoms in 
tetrahedral units. These features were also found in MA-a-GesoSe^. The peaks at around 6.6 and 7.4 A, which are 
also seen in MQ-GeSe2 samples, can be associated, following Ref. [l^, to distances between Se atoms belonging to two 
adjacent tetrahedral units. 

The coordination numbers show in Table [I] were calculated considering the RDF^ MC (r) functions shown in Fig. 0] 
The integrations were made using the following ranges: from 1.7 A to 2.95 A to the first peak and from 2.95 A to 
4.5 A to the second peak, for all RDF^ MC (r) functions. The interatomic distances are also shown in Table [|] For a 
comparison, Fig. EJalso shows the RDF?- (r) functions for MA-a-Ge3oSe7o, and its structural data are also shown in 
Tabled It is interesting to note that except for Ge-Ge pairs, coordination numbers are always larger in MA-a-GeSe4 
than in MA-a-Ge3oSe7o- If we compare our results with those found by MD simulations of ^-GeSe^Si, in which it 
is found A^Qe-Gc = 0.06, A^Ge-Se = 3.87 and A^e-Sc = 1-04, a large increase in the homopolar coordination numbers 
is seen, suggesting that the liquid and the amorphous alloy produced by MA have important differences from the 
structural point of view. This was also verified in an amorphous GasoSeso alloy produced by MA studied recently by 
usSS. 




0123456789 10 



r (A) 

FIG. 4: RDFg^g c (r), RDF§^s e (r) and RDF|" s c e (r) obtained from the RMC simulations for MA-a-GeSe 4 (dashed lines) and 
MA-a-Ge3oSe7c 7 (solid lines). 
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TABLE I: Structural Parameters obtained for MA-a-GeSe.4. 



RMC 







Fir 


it Shell 




Second Shell 


Bond Type 


Gc-Gc 


Ge-Se 


Sc-Gc 


Sc-Sc 




Gc-Gc 


Gc-Sc 


Sc-Gc 


Sc-Sc 


N 


0.22 


3.61 


0.90 


1.71 




2.70 


8.3 


2.1 


10.1 


r (A) 


2.35 


2.36 


2.36 


2.33 




3.80 


3.84 


3.84 


3.75 








MA- 


a-Ge3oSe 7 (^ 










Bond Type 


Gc-Gc 


Ge-Se 


Sc-Gc 


Se-Se 




Gc-Gc 


Ge-Se 


Se-Ge 


Se-Se 


N 


0.26 


3.50 


1.75 


1.25 




3.85 


7.4 


3.7 


9.7 


r (A) 


2.33 


2.35 


2.35 


2.33 




3.83 


3.84 


3.84 


3.75 


MQ-GeSe 2 studied by NO&ii 


Bond Type 


Gc-Gc 


Ge-Se 


Sc-Gc 


Se-Sc 




Gc-Gc 


Gc-Sc° 


Se-Ge a 


Se-Sc 


N 


0.25 


3.7 


1.8 


0.20 




3.2 






9.3 


T (A) 


2.42 


2.36 


2.36 


2.32 




3.57 






3.89 








2-GeSe 2 


studied by ND£ 










Bond Type 


Gc-Gc 


Ge-Sc 


Sc-Gc 


Sc-Sc 




Gc-Gc 


Gc-Sc 


Sc-Gc 


Se-Se 


A r 


0.25 


3.5 


1.7 


0.23 




2.9 


4.0 


2.0 


9.6 


r (A) 


2.33 


2.42 


2.42 


2.30 




3.59 


4.15 


4.15 


3.75 



"These numbers are not given in Refs.0 or0 due to the difficulty in defining the second shell. 



2. Partial Structure Factors 

The partial S% MC (K) are shown in Fig.0 together with the (K) found in Ref. 3 and also with the S% MC {K) for 
MA-a-GezoSeKjL. <Sq^q q (K) has its first three peaks at about 1.1 (very weak), 1.9 and 3.4 A -1 and two minima at 1.2 
and 2.8 A -1 . Their positions agree reasonably well with those found for MQ-GeSe2 studied by ND 3 ^, but intensities 
are very different. The first peak is lower in MA-a-Ge3oSe7o and almost unresolved in MA-a-GeSe4. This peak is 
associated with the FSDP in the S(K) shown in Fig. [5] and, since the FSDP is known to be strongly dependent on the 
Ge-Ge correlations and, to a lesser extent, on the Ge-Se correlationsi^ 7 ^, the low intensity FSDP experimentally 
observed in Fig. [3 could be caused by weak Ge-Ge and Ge-Se correlations at its position. In addition, the heights of 
the second and third peaks are almost the same, and in MQ samples the height of the second peak is twice of that of 
the third peak. 

^Ge^SeC^O nas t w0 maxima at 1.1 and 3.5 A" 1 and a minimum at 1.7 A -1 . Again, the first peak of MA-a-GeSe4 
is lower than that found in MQ-GeSe2 samples, but its position is the same, and it is also associated with the FSDP 
in S(K). The second peak and the first minimum in MQ samples are found at 3.5 A -1 and 2.1 A -1 , respectively, 
indicating that in our alloy they are dislocated to higher (the maxima) and lower- K (the minimum) values. These 
facts can be explained by the important differences between the G™£(r) function and that of the MQ samples for 
r > 4 A. 

Sg^^(K) has two peaks at around 2.0 and 3.6 A~ x , and there is a minimum at 2.8 A" 1 . At about 1.2 A -1 there is 
a small peak associated with the FSDP in S(K). In MQ-GeSe2 samples the peaks are seen at 0.95, 2.05 and 3.55 A" 1 , 
and there is a minimum at 2.75 A~ x . These data indicate that (K) is similar to that of MQ-GeSe2 samples, 

at least in the \ow-K region, concerning peak positions. However, we should note that their heights are different, in 
particular the intensity of the peak at 2.0 A -1 . 

It is important to compare our results considering the BT formalism. Figure shows the S^^ C (K ), S^^ C (K) and 
S^^ C (K) factors obtained using Eqs . ^ O and [31 together with those factors found by ND and shown in Ref. 0- As 
expected, S™ C (K) resembles the XRD S(K) because the scattering lengths of Ge and Se are almost the same, and 
it is very similar to the S^(K) found for MQ-GeSe2, except for the FSDP intensity at 1.0 A -1 . Although there are 
differences in the peak intensities, S^^f c (K) resembles that found for and for MQ-GeSe2— , including 

the sharp minimum at 2.0 A -1 . Scc(K) also behaves like that obtained for £-GeSe 2 ^i and MQ-GeSe 2 -. In Z-GeSe 2 
a sharp FSDP is clearly seen in Scc(K) at around 1.0 A -1 , and this fact also occurs in MQ-GeSe2 samples. MD 
simulations of /-GeSe2 15 ' 16 and MQ-GeSe2 — could not reproduce well this peak. In our case, a very weak FSDP can 
be seen at 1.1 A -1 as a shoulder of the high peak at 2.0 A -1 . Remembering that the FSDP in S(K) (see Fig- EJ nas 
a low intensity and considering all the results discussed above, we should not expect that the FSDP in Scc(K) was 
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FIG. 5: iSGe-Ge(^0, <5Ge-Se(Jf) and iSse-Se(J^) factors obtained from the RMC simulations for MA-a-GeSe4 (dashed lines) and 
MA-a-Ge 3 oSe 70 (thick solid lines), and also the ND Sij(K) of MQ-GeSe 2 (thin solid lines, Ref. [|. 

as high and well defined as it is in I- or MQ-GeSe2 samples. The IRO in MA-a-GeSe4 is different mainly because of 
the introduction of Se-Se first neighbor pairs as "bridges" between the tetrahedral units, thus decreasing the number 
of ES units and increasing the number of CS units. This affects the short and medium range order, which, in its turn, 
changes the concentration-concentration BT factor. 

3. Bond-Angle Distribution Functions 

By defining the partial bond-angle distribution functions Q%ji(cos9) where j is the atom in the corner we calculated 
the angular distribution of the bonds between first neighbor atoms. The six 0y/(cos0) functions are shown in Fig. 
All these functions were calculated considering as r max the position of the first minimum after the peak of the first 
shell (r max « 3.0 A). 

The Oce-Ge-Gefcos 9) function is very noisy because of the very small number of Ge-Ge pairs in the first shell, but it 
shows a tendency for angles around 100°. The 0g c _g c _s e (cos 9) function has peaks at 55-61° and a broad distribution 
from 99 to 118°, with a maximum at 103°. The internal Se-Se-Se angles in perfect tetrahedra are found at 60°, the 
Se-Se-Se angle in trigonal Se is seen at 103° and the angles in small Se chains and rings 46 can be found at 90-116°. 
Thus, the Osc-Sc-Sc(cos 9) function indicates that the tetrahedra in MA-a-GeSe4 are slightly distorted and Se chains 
and rings are formed, in agreement with RS results and with the previous analyses of the Gy (r) functions. The 
@Se-Ge-Se(cos#) function is very similar to the Os c -Sc-Sc(cos 9) function, showing peaks at 58° and 105°, which is close 
to the ideal tetrahedral angle of 109°. 

The 0G e -Se-Ge(cos 9) function peaks at about 58° and 106°. The Ge-Se-Ge sequence in ES units has angles around 
80° that are not seen in this function, reinforcing the small quantity of these units. On the other hand, this sequence 
in CS units has angles around 100°, and n-fold rings have angles ranging from 92° to 125° The Os e -Ge-Gc(cos 9) 
and 6Ge-Se-Se(cos#) functions are similar to the others, having peaks at 58°, associated with threefold rings, and 
at 116° (Se-Ge-Ge) and 109° (Ge-Se-Se), which are related to tetrahedral angles and n-fold rings. The 0jj/(cos0) 
functions above confirm that distorted tetrahedral units are formed in MA-a-GeSe4, with a clear preference for CS 
units. These units seem to be connected by Se-Se bridges, forming small chains and rings, as pointed out by the RS 
data. 
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FIG. 6: Bathia-Thornton Snn(K), Snc(K) and iScc(^) factors obtained from RMC simulations for MA-a-GeSe4 (dashed 
lines) and MA-a-Ge3oSe7o (thick solid lines) and also the ND BT structure factors of MQ-GeSe2 (thin solid lines, Ref. 3). The 
arrow indicates the FSDP in Scc°(K). 




FIG. 7: Qiji(cos0) functions obtained from RMC simulations for MA-a-GeSe4 (squared line) and M A-a-Ge3oSe7o (solid line). 



CONCLUSION 



To summarize, we can conclude that the amorphous GeSe4 alloy can be produced by MA starting from the elemental 
powders of Ge and Se, but the structure of the alloy is different from that found in MQ-, VE- or MG-GeSe2 or MQ- 
GeSe4 samples, making clear the importance of the preparation technique. Structural units similar to distorted CS 
and ES tetrahedra are formed, with a preference for CS tetrahedra, as indicated by Raman spectroscopy. These 
units seem to be connected by Se-Se bridges, as suggested by the high number of these pairs in the first shell, by the 
increase in the distance between Ge-Ge second neighbors, by the @iji(cos 9) functions and also by RS vibrational data. 
These differences in the SRO affect the IRO of the alloy, and this causes the low intensity of the FSDP in S(K) of 
MA-a-GeSe4 when compared to the MQ-GeSe4 alloyi The low intensity of the FSDP in S(K) can be traced back to 
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the partial SG O -Ge(-?0 factor, whose FSDP is related to the ES units, which are found in a small quantity in the alloy. 
Scc(K) reflects these features and shows a very weak FSDP when compared to the factor found for l-GeSe^i^^ or 
MQ-GeSegi. 

As a second remark, this study reinforces the relevance of using and combining RMC simulations with other 
techniques to model amorphous structures, since all features described above were obtained considering directly the 
experimental S(K) in the simulations. 
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